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Various problems regarding the aerodynamics 
of lifting wings are solved. Two methods 
are proposed for replacing the wing, both 
involving "viscous" edge vortices. The 
applications give results which agree well 
with experiments. Two new methods are also 
proposed for calculating propellers based on 
the vortex model consisting of an edge vortex 
and a "viscous" hub vortex. 


12. Kmj (S«1c«ftd trx AaUh3r(i^ 


* 

I 


t*. V.t.ritr Cl.l.ll, (.1 

, Unclejsaif ic4 


IB, OfaiHlwHta 

unclassified and unlimited 


23« Cf«»»4f. (»# 

Tt* Ka. WPas.. 

22. Pmm - 

ITacl&aolfled 

80 



ii 




( 941 rG (tUftL- 


-) 


A/ -/ 5 - 3 ^ 9 3.1 
A/r3- 315 - 7 ^ 
K! -Ka 91 / 






CONTRIBUTION TO THE AERODYNAMIC STUDY OF WINGS AND PROPELLERS* 
Maurice Menard and M. A. Toussaint 

PREFACE 

In order to develop an experimental method to study the aero- 
dynamic functioning of profiles of propeller blades, M. Menard 
first solved various problems about the aerodynamics of lifting 
wings. 

First of all, he proposes A calculation method for the correct- 
ions caused by lateral boundary layers for "between panels" experi- 
ments, using the logarithmic law for the velocity distribution. In 
boundary layers and for the correlation variation of the circulation 
over a section of the experimental wing, he calculates the lodal 
induced angles and the average Induced angles for this wing section. 
In this way, we obtain the corrections to be applied to the unit 
lift curves in order to extrapolate them for infinite length. 

Various applications perfectly justify the use of these corrections 
which occur for "between panel" tests. After this, M. Menard makes 
an Important and original contribution to the calculation of finite 
span wings. He proposes two methods of replacing the wing which both 
Involve "viscous" edge vortices. This new model seems to represent 
real edge vortices very well. The applications made for the median 
section of the wings, and those for any other sections, give results 
which agree well with experiments, especially as far as the "edge 
effect" is concerned for the distribution of the lift force over the 
span of rectangular wings. 


m 
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For the experimental study of the aerodynamic functioning of 
the profiles of a propeller blade, M. Menard considered and devel- 
oped a method of measuring the local aercrdynamlc Incidence. This 
method consists of measuring the aerodynamic pressures at 3 points 
located close to the leading edge of the blade profile. The trans- 
mission of these pressures to micro-manometers required a "celling 
rotating Joint" having 3 compartments. 

Once the aerodynamic incidence of these various profiles of a 

blade are. known and the characteristics of the same profiles for 

infinite span are known, one can calculate the radial distribution 

of the traction and coupling coefficients. By Integrating this,, one 

• V 

obtains the global coefficients for various values of Y = ^ • The 

results calculated In this way agree well with the global measurements 
made with a balance. 

Finally, M. Menard proposed two new methods for calculating 
propellers based on the vortex model consisting of an edge vortex 
and a "viscous" hub vortex. 

The application of these methods to the particular case of a 
propeller sr.ov;s good agreement between the calculated and measured 
Induced velocities. 

In conclusion. It Is obvious that the work of M. Menard will 
Interest aerodynamlclsts In laboratories and in study offices. In 
closing, I would like to confirm that I appreciated the collaboration 
at the Salnt-Cyr Aero ’ Technical Institute. 


A. Toussalnt, 

Professor at the Scientific 
Faculty, Director of the 
Aero Technical Institute 
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Chapter I 


A* 

GENERAL REMARKS 

Modern calculation methods for propellers use the principles 
described by Joukowski [1] in 1912. 

At this time, the calculation methods were based on the applica- 
tion of theorems of momentum and the moment of the momentum. These 
are theorems which allow one to define the values of the force and 
the. couple of the propeller. 

Joukowski attributes additional axial, radial and rotational 
velocities near the propeller to vortices, the existence of which 
was demonstrated by Plamm. 

I . Propeller with arbitrary circulation distribution 

In this case, the model of the vortex system used by Joukowski 
consists of the following: 

a) radial vortices tied to the blade located in the propeller 
rotation plane; 

b) free vortices which depart from the blade and rollup into 
solenoids. 

The circulation of the vortex tied to the blade is defined as 
a function of the local lift coefficient, the relative resultant 
velocity to which the profile is subjected, and the length of the 
profile chord at the radius under consideration: 



r, circulation; 

C^, local lift coefficient; 

• W, resulting velocity at the radius under consideration; 

* 

Numbers in margin indicate pagination of foreign text 


5 



1, chord of the profile. 

Motion Is assumed to take place In a perfect fluid, and the 
drag coefficient does not appear In this theory. 

II. Propeller with constant circulation distribution /2 


Joukowskl gives a constant value of to the bound vortex. 

The vortex sheet which makes up the wake of the propeller Is replaced 
by an axial vortex having circulation or 2 F^ (two-blade pro- 
peller) and a vortex at the tip of the blade with circulation F^ 
which rolls up Into a solenoid. 

The velocity Induced by an element of the vortex of given 
length Is defined by the Blot and Savart electromagnetic Induction 
law. 


For propellers with arbitrary circulation distribution, 
Joukowskl points out the difficulty Involved In calculating the 
values of Induced velocity for any radius of a propeller using 
three hellcoldal vortices. 

In order to avoid this difficulty, the number of blades Is 
Increased Indefinitely. This allows one to define average Induced 
velocities for any radius of the propeller blades. 

This calculation of the Induced velocities allows one to define 
the angle of attack of the propeller for any section under consi- 
deration. This section Is supposed to have the characteristics of 
an Infinite span wing (same chord) calculated for the same aero- 
dynamic Incidence. 

By Indefinitely Increasing the number of propeller blades, 
Joukowskl substitutes a fictitious flow for the real flow which has 
on the average the properties which are comparable with those of 
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real flow. 


III. Extension of the Joukowskl theories 


The theories derived from the Joukowski theory can be separ- 
ated into two categories: 

a) theories which use an infinite number of blades. 

This involves a continuous vortex distribution. The induced 
velocities at radius r have a constant value along a circumference 
of radius r (average induced velocities). 

A correction has to be applied in order to take into account 
the finite number of blades of the real propeller. 

b) theories which use a finite number of blades. 

This is proceeded by a calculation of the velocities induced 
at a given point of the blade by vortex strips which emanate from 
the blade and roll up into solenoids. The Induced velocities have, 
therefore j been determined and are called "local” velocities. 

/3 


PURPOSE OF THE PRESENT STUDY 

Numerous calculation methods for propellers have been proposed. 
The comparison with experiments shows that for propellers operating 
at high similitude — ^ values, these methods provide sufficiently 
accurate approximations . It seems that there is a large disagree- 
ment between theory and experiment when the propellers operate at 
small values of v = — 

The purpose of this study is to research measurement ^procedur^ 
which will allow the determination of experimental conditions of 
aerodynamic operation of profiles which are used for propeller blades 
with a good degree of approximation. For this purpose, two 
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dT dR. ^ dR, . 

_ = _cos9- — sm9 


— P ^ W" / (C; COS 9 — Cj sin 9), 

dQ /dRr I dR- . \ 

= P 5 W* I (Cr COS 9 + C; sin 9), 



measurement methods were used: 

- measurement of an average velocity 

- direct measurements on the blade. 


I . Measurement of average velocity 

This is an an anemometric study of the blowing of a propeller, 
a study which was performed with the usual measurement procedures 
(directional probe). 

II . Direct measurements on the blade 

Let us consider a propeller blade profile located at radius r; 
the resultant of the forces applied to this profile Is equal to the 
geometric sum of lift and drag. 

Let W be the resulting velocity over the blade at the radius 
under consideration. We can write: 

dR^ = I C.W*/ dr. 

dR^ = He, W*/ dr. 

In Figure 1, It Is possible to describe the traction and the 
elementary couple for radius r. 

/i 

The relationship between the angle cj) and the geometric step 
angle 9 of the propeller Is: 



(}) = e - i; 

1 = Incidence of the profile under consideration. 

The angle 4>, Incidence 1 and resulting velocity W determine 
the operational conditions of the propeller. 

MEASUREMENT OP THE INCIDENCE AND THE VELOCITY 

Let us consider a wing profile placed In a flat stream In a 
wind tunnel. Over this profile, we measured the pressure at three 
points. The pressure differences between two of these points and 
the third point (Ap^ and AP2) are proportional to the square velo- 
city of the flow at Infinity and at the aerodynamic Incidence angle 
of the profile. 

Since we have a system of equations with two unknowns, It Is 
possible to determine the aerodynamic Incidence of the profile and 
the velocity at Infinity, under the condition that we have first 
made pressure tap measurements. We use this procedure In order to 
determine the aerodynamic Incidence and the velocity resulting at 
various radii of a propeller blade. 

In order to measure these pressure taps, we constructed a wind 
tunnel with a flat rectangular test section In which we wish to 
directly measure the pressure taps over the blade. But the varia- 
tion of the chord and the twist of the blade, for a section sub- 
jected to the flow, made these measurements Inaccurate (variation of 
the circulation along the stand). 

Under these conditions, we then studied pressure taps over a 
cylindrical wing with a constant chord and which had no twist. In 
order to construct a propeller having the same base profile. The 
lift of this wing was measured by Integration of the pressure along 
the chord (median section) . The slope of the unit curve of lift 
determined under these conditions was too small. This showed that 
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the fiotltious equivalent aspect ratio obtained for plane flow (flat 
wind tunnel) did not have the characteristics of an experimental 
infinite aspect ratio, obtained between panels and with a geo- 
metric aspect ratio of 8.66 (special measurements performed in the 
median section of the wing) . 

V/e tried to establish the causes of this anomaly. Prom this, 
we then studied the influence of boundary layers which developed 
along wind tunnel walls for plane flow tests. 

It seemed to us to be necessary for possible application to 
propeller blade measurements to know the Influence of the aerodynamic 
flow field curvature close to a wing having finite span (curvature 
due “o the presence of Induced velocities which vary along the chord 
of the profile). 

In the thesis of M. Carafoll, we note the difficulty encountered 
by the author when trying to superimpose the diagrams of pressure 
measured in the median section of wings having varying aspect ratio 
(diagrams reduced to the same aerodynamic incidence). This document 
led us to use a new model for the vortex system of a wing with finite 
span. 

/5 

The new replacement system proposed allowed us to do the follow- 
ing: 

a) calculate the distribution of the circulation, along the 
span of a wing having finite aspect ratio; 

b) to obtain the superposition of the diagrams for the median 
section of different finite aspect ratio wings, taking into account 
the field curvature. This was also done for wings tested between 
panels and which had different geometric aspect ratios; 

c) to correct measured pressures measured over the propeller 
blade (in order to take into account the fl^d curvature); 

d) to establish a new vortex model whlcri'ds simple and which 
can be used to calculate a propeller. 
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Chapter II 


STUDY OP AN EXPERIMENTAL INFINITE ASPECT RATIO 


I. Definition 


In aerodynamics, one usually considers a cylindrical wing which 
is displaced perpendicular to its generators in a perfect fluid. 
Under these conditions, one has realized a theoretically Infinite 
aspect ratio. The motion of the fluid around the wing is called 
"two-dimensional" or "plane flow". The flow is identical for all 
sections perpendicular to the wing generators. 

It is very difficult to realize a wing whose geometric aspect 
ratio is very high (pseudo-infinite). All of the experimental con- 
figurations built to obtain a correct experimental infinite aspect 
ratio attempt to make the results independent of the geometric 
aspect ratio of the wing under study. 

Among all of these realizations, we will now consider the one 
used at the Salnt-Cyr Air Technical Institute (test between panels). 

Professor Toussalnt [2] has stated that infinite experimental 
aspect ratio is difficult to establish in the wind tunnel for the 
following reason: 

- The transverse distribution of velocity between the two par- 
allel walls of a wind tunnel is not flat. Along these walls, bound- 
ary layers develop. The motion of the fluid particles around the 
profile under consideration is not uniform along the generators of 
the wing under study. Therefore, we no longer have a "plane" flow. 

II . Theoretical study of the ppeudo-plane flow 

Let us consider a wing A of span L having a constant profile 
along the span and not warped placed between two parallel walls PP* 



cf a wind tunnel through which a fluid current flows of velocity V. 

We will study the influence of boundary layers existing on the walls 
PP’ on the aerodynamic characteristics of the wing (Figure 2). 

For edge regions of the wing located inside the boundary layers 
B, C, the velocity of the flow is no longer V but it is U. 

/8 

The pseudo-plane experimental flow will be considered as being 
equivalent to the linear juxtaposition of the following: 

a) a theoretical plane flow for which the aerodynamic incidences 
are the same over the the entire span of the wing (geometric inci- 
dence referred to the direction of 0 width of the profile) 

b) a flow caused by the velocities induced by the system of 
free vortices coming from edge regions where the circulation is 
variable due to the fact that the velocity U is different from V. 

The juxtaposition of these two flows leads to a system of the 
integral-differential type for determining the velocities induced, 
along the span of the wing. In order to simplify the calculations, 
we assume that the variation of the circulation along the span of 
the wing is only caused by the transverse gradient of the velocity 
which exists in the boundary layer. This amounts to ignoring the 
action of flow b). 

Let : 


i, geometric incidence in the median plane; 
dc 

z, slope of the unit curve of lift for the median section; 

dl 

1, chord of the wing. 


The circulation in the median section has the following value: 


( 1 ) 





di 2 
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The circulation In the plane having ordinate value y Is equal to: 


1 . Diagram of the equivalent vortex system 

Just like in the case of a wing whose span is ^'inite (therem 
of Prandtl), we assume a single vortex placed at the focus of the 
wing, a vortex which has the same circulation as the real wing. 

/9 


The attached vortex has the same transverse dimension as the 
wing and the same circulations as the real wing at the- same points 
at each of the points. 

Let us consider the extremity of the wing over which a circula- 
tion gradient Is established (Figure 3 ) • 

In the plane with the ordinate y, the circulation Is f . In the 

y 

plane y + dy, the circulation Is: 


For the planes with ordinates y and y + dy, therefore, we have 
the circulation: 

This elementary circulation will be considered as the circula- 
tion of the free vortex emanating from the profile at this point. 

The collection of the free vortices will constitute a plane 
sheet . 

Also, the flow exists between walls and It is necessary to take 
Into account the stresses Imparted to the flow by the walls. 
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Fig. 3 

In the replacement vortex system, we represent these stresses 
by Infinite rows of vortices with alternate orientations. 


If we consider the free vortex with circulation emanating 

from the profile In the plane y^ the vortex system which takes Into 
account the presence of the walls Is the one described by Figure 4. 

By symmetry. It Is obvious that at a point with an ordinate 
-y, a free vortex with the circulation +df will detach. Taking 

y 

Intc account this remark, the diagram of the vortex system becomes 
the one shown In Figure 5. 


We are dealing with an Infinite row of vortices with alternate 
orientations ~dV,j of step L and with ordinate origins of ±y . 


The complex potential for the row +df Is: 

y 


/j (r) = _ t ^ S Log (r -r y) + 2 ^ 


For the row df , we have 

y’ 


fz ^ 1 Log (; - y) -f 51 Log [(: - y)' — m- L-]. 


For the two rows, the complex potential Is then: 




Log 


(: 




If) \(z -r y)- -f L^1 f (: + y-) + 2 U-T f ■ ■ .] \(z + y)» + m - LM 
y) [(-’ - y)^ - I {- - y^) - 2 L^] - y)’- - L^r 



By making an expansion, we obtain the formula of M. Vlllat : 


p, — 

- * . 77 , . 

sin-(r — y) 

The complex velocity corresponding to this will be: 


dw 


= ‘ ^ I £ (- -f- y) — cotg £ (i — y)] 




.dr,, 


sin 


2jl2 


L 2nw 2-; 
cos — — cos —|— 


At a point with an ordinate Y, the complex velocity is reduced 
to the vertical component and its value Is given by: 


Al 


(4) 


dv — 


dTy 

L 2 :t 1/ 

CO® L 


• 2 71 1/ 

sin — 


COS 


27 tY 


Infi 
(er.a 
to t 

(5) 


In the calculation we have assumed that the free vortices were 
r.ite in the two directions. But here we only have semi-vortices 
natlng from the wing) and by symmetry, the expression (4) has 
e divided by 2: 


. dT„ 



cos 


2jLy 

L 


— cos 


2 t; Y * 
L 


The induced velocity resulting from the point Y will be equal 
to the sum of the velocities Induced by all of the free vortices and 
their Images at, point Y, that is: 


L 



2 . Representation of velocities in the boundary layer 
a ) Friction law 
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The formula generally used for friction of a flow over a flat 
plate without a longitudinal gradient of static pressure is the one 
given by Prandtl: 


( 6 ) 



with 



Let us apply the theorem of momentum to the boundary layer. We 

have : 

( 7 ) ^ ^ 

.pVo* dx' 


This expression allows one to express the friction as a function 
of Reynolds number. If we consider the equations (6) and ("), we 


or ; 


(fS“ 

dx 


0.01255 


1’ 



dp 

dx 


0.01255 

1 
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Let us Integrate this expression between the limits 0 and x. 
We have: 

4 _ 0.01255 X 

(W 

from which we find 


( 8 ) 


0,036 

X l’ 
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with 


and finally 




0,0288. 
1 » 

(yl.T)5 


b ) distribution of the velocities In the boundary layer 


In the boundary layers which developed over a flat plate sub- 
jected to a turbulent flow without longitudinal static pressure 
gradient, the velocity profiles are similar, that Is, they can be 
represented by a single curve. 

We can represent the velocity profile In the boundary layer 
by using a lav; of the ' following form: 


( 9 ) 


i-l+K,Logg). 


( 10 ) 


ror the following calculations, this leads us to the following: 

.j’'KLog(2)„(S) = K,*. 

” s) (0 (s)’ 


( 11 ) 


•s« = K, (1 — 2 Ki) S. 



/ 1 3 . 


If we consider equations (10) and (11), we can see that It Is 
possible to obtain by graphical integration the quantities 5^ and 
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XX 

5 If we have an experimental distribution of velocity in the bound- 
ary layer. This allows us to determine the value of the coefficient 
introduced in formula (9). M. Michel [3] shows that with this 
coefficient it is possible to determine a thickness 5^ which is not 
the ordinate where u = V but which determines in a more accurate 

X XX 

manner the values of 6 and 6 by associating them with the values 
of K^. The friction values and the values of flow rate are therefore 
represented better. 


Experimentally, we have noted a small deviation between the 
thicknesses 6^^ and 6. 

We can relate the shape parameter H to the coefficient K^: 


1 - 2 k; 

Our representation of velocity (9) can be written as: 



H — 1 
2 H 



For a flat plate without a longitudinal static pressure grad- 
ient in turbulent flow, the value of the shape parameter H is equal 
to 1.3 j that is: 

Ki = 0,116. 


We finally find the following as a representation law of velo- 
city in the boundary layer: 

(12) 1+0.116 Log (a). 

li * y 

3. Introduction of the law ^ = l + KiLog| as a representation of 
the the boundary layer * • . 

It is useful to define the transverse circulation gradient in 
order to find the velocity Induced at an arbitrary section of the 
wing . 
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We have seen that It can be written in the form: 




Formulas (1) and (2) of paragraph II result in the following: 


(13) 







u 

Vo' 


In the preceding paragraph we have seen that it was possible 
to adapt formula (9) as a velocity law in the boundary layer from 
which we find 


/in 


(14) 


r; = ro(i + KiLog|j. 


By considering equations (8) and (11), we find 


1 

Ki (I - 2 Ki)’ 


— = 11,22 


with Ki = 0,116, 


(15) s - 11.22 X 

This latter formula allows one to calculate the thickness of 
the boundary layer at the measurement point x (this is only true in 
the case where we assume turbulent flow over a flat smooth plate 
without any static pressure gradient in the longitudinal direction) . 


0.036 

\0,2 


Let us consider figure 6; in the formula (l4), ^ will be 

defined by: 

- L 

? = 2 " y- 

The equation (14) can be written 
as : 


(16) 


T: = r, / I + Kx Log I - 




The transverere- circulation gra- 
dient is: 

(17) dr, = ^dy = Kir,— ^ 





'S 
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We have found in paragraph 1 that the velocity induced at a 
point with ordinate Y by a free vortex and its images was: 


( 18 ) 


do = 


sin^ 


dr^ 

2 L 2 X w 2^Y' 

cos -r-^ — cos 

L# L* 

By relating equations (1?) and (18), we have 


(19) 


du==^- 

2L 


sm 


Ijlh 


Ail 


(«- t ) 

The resultant velocity Induced at the point with ordinate Y 
by vortices ’’from the boundary layer” and their Images will be 
given by: 

ri ■ 

(20) " 2L 




Ail 


(cos cos (y-|) 

2 

and the integration can be performed graphically. 
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The velocity diagram at the ordinate point will be the one 
shown in Figure 7- We can rewrite: 


^ = tgSf^^Siv 

' It 





Figure 7 


V (axial velocity at distance Y 
«y 

from the axis). 

Two cases have to be considered: 


a) the ordinate Y is located inside the boundary layer, the 

velocity of the flow is equal to V and the local induced angle can 

be written as given by (20) and by using ro = c,^ , we have: 

2 


( 21 ) 


Sir. 


K. C,l 
4L 
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b) the ordinate Y is located inside the boundary layer, the 
velocity of the flow at this point a is U, and the local induced 
angle becomes: 


( 22 ) 



It is possible in equation (22) to use the logarithmic represen- 
tation of the velocity in the boundary layer: 

K, C,l 


(23) 


S j'u = 





JlL 




We will define an average Induced angle as follows: 
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(24) 



III . Calculation of the aerodynamic characteristics of a wing 
placed in a pseudo-plane flow 


We have defined an angle Induced at the section Y of a wing 
placed in a pseudo-plane flow. Because of this Induced angle, the 
usual reference axis system (directions parallel and perpendicular 
to the velocity direction at Infinity upstream) will be changed. 

In Figure 8, we see that the resultant of the forces applied 
is no longer perpendicular to V but perpendicular to V. We then 

O' 

have : 
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dRj = dR cos S iy 

= pVyTydy cos 8 iy 
= pVj,r„dy^pVrj,dy. 

dR.^dR. 


Thus, the modification of the lift can be considered as small. 
On the other hand, as far as the drag Is concerned, we show a very 
substantial Increase which can be written as: 


since : 


dR* = dR sin 8 ly = ? Ty Vy sin 8 iy dy 

= ?TyVydy'^^dR: 

dR* = ^ dR„ 


(25; 



dR* = 8 iy dR,. 


We can determine the Induced local resistance coefficient: 


(2f 

(27: 


dR, = C,^V*ydy/; 

dR* = C*^ V*ydy/; 


using (26), (27) and (25), we have: 

C*. . I V*y / dy = 8 fyC, I V*y/dy. 

Cj., = C, 8 iy. 

/1 7 

In paragraph II, we have defined an average Induced angle from 
an examination of the distribution of induced local angles along 
the span of the wing. The Induced angle had the following value: 



The lift and drag coefficients measured with aerodynamic bal- 
ances are average coefficients (since we weigh the entire wing). • 

In Figure 9, we show the unit curves of lift. I represents the unit 
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curve of lift measured with the balance. 


We have seen that for a given 
valui of the lift coefficient, the 
reference axis position changes. 


For a lift coefficient value of 

, the aerodynamic Incidence 
Is therefore not a but a - Al. In ^ 

order to obtain the corrected unit 



I 




I 

1 

r 

1 

1 


curve of lift. It Is sufficient to 



1 

J 

1 

(• 

plot the following curve: 

C, = /(* — Ai). 

Ck&t 

t'ig. U 

CC 


Before we have seen that there was an Increase 
drag. The drag coefficient measured with a balance 

In 

can 

the 

be 

"Induced 

written 


as : 

Cii, = Cx, + Cx.„," 


The Indices: 

B = balance 
C = corrected 
1 = Induced 
have been defined. 

The Is equal to: 

c -c'-c 

In order to make corrections to 
the unit drag curves, C = F(l) 
measured with the balance, one has 
to take Into account the fact that 
there Is an Induced drag and that 
the geometric Incidence of the wing 
Is reduced by an amount A1 (average 
Induced angle). 





In Figure 10, we show successive corrections to be applied to 
unit curves of drag measured with the balance. The distance AB 
represents the Induced for a value of C_ corresponding to the 
incidence a. In order to take into account^the average Induced 
angle, it is necessary to attribute a new value C derived from (30) 
to the aerodynamic incidence (x - Ai). The distance BC therefore 
represents (x - Ai). 
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Remark . When we measured the drag coefficient by the wake 
method, we only measured the value Cx? . in order to make correct- 
ions, it is. necessary to take into account the local Incidence 
variation due to the local Induced angle (measurement plane). 

IV. Experimental configuration 

1 . Determination of lateral boundary layer characteristics 

Before each series of tests, we made velocity measurements In 
the lateral boundary layers of the wind tunnel. We made these mea- 
surements essentially at the location where the focus of the wing 
was located. In order to perform these measurements, we used a 
total pressure tube whose extremity was flattened and the sides 
were defined as follows (Interior): 0.10 mm x 2.5 mm; exterior: 

0.5 mm X 3*5 mm). This total pressure curve was closed off by a 
carriage which allowed displacements of 0.2 mm. 

2 . Measurement of aerodynamic forces (lift) 

We used three measurement procedures: 

a) Brass tubes Inserted in the profile. 

The profiles were made of walnut. In the profiles we engraved 
, grooves parallel to tl^ generators using a mill. Brass tub”esj 0.7 x 
1.5 mm, were Introduced in these grooves and were glued using strong 
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glue. We provided for a correct surface of the profile by using 
cellulose mastic In order to cover up the roughnesses caused by 
milling of the profile. 

The pressure probe holes had a diameter of 0.5 mm and were 
arranged over the same section of the wing In order to avoid measure- 
ment errors Introduced by a fast variation of the Induced angle along 
the span. 

b) Pressure measurements on the profile using a licking probe. 

For this purpose, we used a probe having special characteristics. 
We did not wish to use the carriage which supports a licking probe 
because this procedure has a certain number of disadvantages. For 
example, It Is luminous and there Is risk of Introducing local per- 
turbations In the circulation around the profile being measured. 

This makes It difficult to measure In regions around the leading 
edge and around the- trailing edge of the wing. 

The measurement probe is a modern type as far as its shape of 
its extremities are concerned. In Figure 11, we show a diagram 
shov;lng the shape of the probe support: 

- AB, pressure probe which can pivot around axis B; 

- BC, small rod which can turn around the axis C; 

- at B and C, reference springs which make the part AB be par- 
allel to the profile and therefore parallel to the flow. 

/19 

During the calibration of this probe, we noted that there was 
an Interaction between the part BC and the part AB. When the distance 
BC Is greater than 20 times the diameter of the stem DC, Indica- 
tions of the probe are Independent of the position of BC. 

With this measurement device, it was possible to measure press- 
ures over wing profiles from the leading edge up to 20 mm from the 
trailing edge. 
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The part DC was connected to a device which allowed the trans- 
lation of the probe parallel to the chord of the profile studied. 

The calibration curve of the probe is shown in plate 1. 


The raeasurement of the pressures over the wing allowed us to 
determine the lift of this wing. With consideration of Figure 12, 
the lift of the wing can be represented by the integral: 

P = / pds cos ot 
J c 

■'* * 

v;hi3h is an integral which extends over the entire contour of the 
prcflle. In order to carry out the integration, we construct the 
pressure diagram as a function of the abscissa OX of the projection 
of she measurement point on the flow direction with a flow speed of 


c) Measurement of the lift using an aerodynamic balance 

In wind tunnel no. 1 at Salnt-Cyr, there is a pseudo-rectangular 
test section 'with a width of 0.78 meters and a height of 2 meters. 

The wing was placed between the panels and is connected by means of 
a wire network to a typ,e* Denls-Gruson aerody^namlc. balance. The 
lateral play allowed between»*t.he wing ancf the. panel is 2 mm. 
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3 . Measurement of aerodynamic forces (drap;) 


We determined the drag of the wings tested either by using an 
aerodynamic balance or by. the wake method. In the latter case, we 
started with measuring the velocity in the wakes of the wings In 
order to find the loss In momentum of the fluid close to the pro- 
files under consideration. We used formulas established by Jones 
which allow one to determine the drag coefficient from static press- 
ure measurements and total pressure measurements performed upstream 
and downstream of the wing: /2Q 

total pressure at Infinity upstream; 
total pressure in the measurement plane; 

Vq, velocity at Infinity upsream: 

AH = H-, - 

1 o 

Vie have; 

/ 

0 

0 and y are wake limits; 

K Is the correction coefficient of Jones, taking Into account 
the differences In static pressure between Infinity upstream 
and the measurement point. 

V. Presentation of results 


In the experiments we allowed the various quantities which 
Intervene in the correction formulas of the preceding paragraphs to 
vary : 

I 

1 . Variation of the characteristics of lateral boundary layers 

We constructed a wind tunnel with a rectangular cross-section 
(0.050m X 0.50m). The test section length is 1 meter (see Figure 13 
and photograph Figure l4). 


28 




We made the longitudinal static pressure gradient 0 by diver- 
gence of the small dimension walls. On the large side of the wind 
tunnelj we drilled orifices which allowed us to Introduce side 
plates carrying the impression of the profile under study. 

The static pressures were measured in the median section of a 
v;lng with a chord of 132 mm (Figure 15) for different incidence 
angles. This allowed us to determine the corresponding C„ values. 

Zi 



Figure 14 


The characteristics of the lateral boundary layers were changed 
with the surface state of the collector. 

The curves (plate 2) show the velocity distributions In the'- 
boundary layers at the measurement point. The wind tunnel collector 
was rough or smooth. 
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Figure 15 


The lateral boundary layers have the following characteristics: 

a) rough collector: 

= 16.5 mm, = 0.115; 

b) smooth collector: 

5t = 12.8 mm, = 0.128 

Lj X 

The curves C„ = F(l) for the two configurations are shown In 
plate 3 and the corresponding values are given In Tables 1 and 2: 

TABLE 1. Rough collector 


i ' 



■ 

2® 

4® 

6® 

8® 

0 

0 

0 

0 

n 

- 27,5 

— 11,1 

T- 18, Q 

-i- 34,4 

o0j3 

-t 65,3 

-f- 78,7 


TABLE 2. Smooth collector /22 


1® 

— 3®, 2 

— i^.io ' 

4®,8 

6®,8 

1 

8^8 \ 

100 C: 

— 21 ,5 

-4,7 ' 

I 

1 

42,5 

60 

1 

74,7 ; 

1 
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The Induced local angle for the median section is determined 

from: 



for a lift coefficient value taken equal to 100 C_ = 50. 

The calculated Induced angles are the following: 

- rough collector: 2.59° 

- smooth collector: 2.07° 

The slopes of the unit curves C = f(i) corresponding to this 

are : 

- rough collector: 0.0757; 

- smooth collector: O.O 815 

The slopes of the curve C = f(l) after corrections are equal 

to : 

- rough collector: 0.1245; 

- smooth collector: 0.1233 

In plate 3 we show the unit curves. We can see that between 
the curves obtained after corrections for the boundary layer effect 
and the other curve, there Is a slight deviation which is on the 
order of magnitude of the experimental errors. 

2 . Variations of the geometric aspect ratio of the profile 

For these test, we used a Clark Y type profile with a relative 
thickness of 18^ and a cord of 1 = 0.090 m. We measure the pressures 
in the median section of the wing profile for two very different 
geometric aspect ratios. The wing is placed in a flat wind tunnel 
with a section of 0 .050m x 0.50m and then in a Lelarge wind tunnel 
with plane flow having a section of 0.78m x 1.50m. In the two cases 
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Figure l6 


considered, the wing traverses the wind tunnel from one wall to the 
other without lateral play (Figure 16). 

The geometric aspect ratios of the wing for two test '.conditions 
are the following: 

- 0.050m X 0.50m wind tunnel 

= 0.55; 

- 0 . 78 m X 1 . 50 m wind tunnel 

X = 8.66 

The height of the plane flow is different in the two cases and 
it is necessary to make a correction to take into account the limited 
height of the fluid channel. We will use an experimental correction 
of M. Tous saint: 

= 1 

f. • P' 

1 -r 0,863 

hr 

h, height of the fluid channel; 

1, chord of the profile; 

Ci— • value of the lift coefficient for an infinite height of 
the fluid channel; 

value of the lift coefficient for a height h of the fluid 
channel. 

The lift coefficients derived from pressure measurements for 
the Clark Y profile are given in Tables 3 and ^ and the corresponding 
C„ = F(l) curves are shown in plate 4: 



TABLE 3. — = 0.55 


1 1 .. 

— 8 

— 4 ° 

— 2“ ! 

0® 

■| 

100 C, 

4 

36 

• 53,5 I 

1 

69 

85,5 , 


TABLE -4. — Xg = 8,66 



1 

— 8® 1 

— 6® 

— 4® 

-2® 

-0®,o0| 

1 

4 

6® 

8 

— .j 
1 

10® i 

100 c. 

7,5 1 

I 

1 

29 

49 

68,5 

i 

81,5 

1 

99 

122 

135 

! 145,5 

1 

158.7 


The characteristics of the. boundary layers for the two tests 
are the following: 

- 0.050m X 0.50m wind tunnel: 


= 12.8 mm 

- 0.78m X 1.50 m wind tunnel: 
0^ = 50.8 mm 


= 0.128 


= 0.12 


The slopes of the unit curves C = F(l) are; 
- 0 . 050 m X 0.50 wind tunnel: 
dC„ 

= 0.0807 


/2 4 


dl 


- 0.78m X 1 . 50 m wind tunnel: 
dC_ 


di 


= 0.103 


The corrections to be made to the above coefficients are the 
following in order to take into account the limited height of the 
fluid channel: 

- 0.050m X 0 . 50 m wind tunnel: 


C; 1 

Cr. “ 1,028 

- 0.78m X 1 . 50 m wind tunnel: 


0 , 973 ; 


C-„ 1 , 0031 * 
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In the latter case, the correction can be ignored. 


For a value of the coefficient 100 C_ = 50, the angle Induced 
in the median section of the wing is equal to: 

0.050m X 0.50m wind tunnel: 

S t„_o = 1°.4I ; 

0 . 78 m X 1.50m wind tunnel: 

S ttf.o = 00,019. 


lay 


The slopes of the unit curves corrected for the 
er effect" are: 

0 . 050 m X 0.50 wind tunnel: 


"boundary 


di 


0.1017; 


0 . 78 m X 1.50m wind tunnel: 


= 0.103. 
di 

We can see on the one hand that we have practically realized 
an infinite experimental aspect ratio in the Lelarge wind tunnel. 

The angle induced for the median section is very small. Also, the 
unit curves for lifts obtained after correction have slopes v/hich 
are independent of the geometric aspect ratio of the wing. 

3 . Determination of the global forces applied to the wing 

In order to perform these tests, we used wind tunnel no. at Salnt- 
Cyr. The profile study had the following characteristics: 

- profile chord 1 = 0.600m; relative thickness 12%, 

‘biconvex; geometric aspect ratio = 1.3. 

The unit lift curve of this profile was determined by measuring 
pressures in the median section. The unit drag curve was determined 
by measuring the wakes (median section). The global forces -which 
apply to the wing were measured using an aerodynamic balance. 



a) Lift measurements 


The curve which gives the velocity profile In the lateral bound- 
ary layers Is represented In plate 5- The determination of the char- 
acteristics of the lateral boundary layers gives: 

= 31.5 mm = 0.1125. 

Using graphic Integration, It Is possible to determine the 
variation of the Induced angle along the semi-span of the wing using 
formula (23). Plate 6 shows this variation for a value of 100 C_ = 

’Zi 

100 . 


The average value of the Induced angle corresponding to this Is 
equal to 


A1 = 1.15“ 


The calculated Induced angle for the median section Is equal to: 

In Table 5j we give the values of the coefficients C obtained 
by Integration of the pressures measured In the median section of the 
wing for various Incidences; 


TABLE 5 


1° 

Qo 

2® 

i 

i ■*° 

6» ; 

8” 

10'’ 

100 c. 1 

1 

15,5 

• 3 R 

06 

1 57,3 

i 

78,7 i 

99 

114,5 


In Table 6, we show the lift coefficients and drag coefficients 
measured with the balance; 

TABLE 6 


% 

1 — 3°,60 

— 1»,60 

0<’,50 

20,50 

40,50 

60,60 

80,60 

100,60 

100 C: 

! - 19,1 

— 2,65 

19,1 

59 

57,9 

77 

os 

112.2 

100 Cx 

i 1,21 

J 

1,05 

1,055 

1,43 

1 

2,02 

2,67 

i 

3,68 

5 !i 

li 
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If we consider plate 6, we can see a difference between the 
slopes of the lift curve determined from pressure measurements on 
the profile and those derived from balance measurements. 

For 100 C = 100, the angular variation of the reference axis 
system corresponding to the slope changes gives an Induced angle 
equal to 1®. It should be noted that the unit lift curve measured 
from pressure measurements In the median section is not the unit 
curve corresponding to an Infinite experimental aspect ratio . The 
calculation shows that there is a local induced angle for the section 
under consideration which Is equal to: 

61 = 0.15° 

y=o 

In order to correct the slope of the lift curve determined with 

* 

the balance. It is necessary to take into account the change In the 
angular position of the reference axis system equal to 

A1 = 1° + 61 = 1.15° 

jr W 

This angle is exactly equal to the average Induced angle determined 
by calculation. 

Therefore, the curve C„ = f(l). resulting from global measurements 
from the balance amounts to the curve C = f(i) for infinite exper- 
Imental conditions. 


Remark . We can determine the equivalent fictitious aspect ratio 
obtained for tests between panels in wind tunnel no. 1 for a wing 
with a chord of 0.600 m, by assimilating the wing between the panels 
into a wing having a finite span with aspect ratio such that: 

= % X 57,3; 

TtX 

for 100 = 100, we have: 


10 1- 57,3 

A I = l«,lo = — — . 
- 

- 3.61 “ ’ 
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b) Drag measurements 


In Table 7, we show the values of measured by the wake 
method for the median section of the wing. 


TABLE 7 


‘1 ’ 
■ • 

— 2® 

0® 

2° 

4® 

6® 1 

8® 

!] 100 Cx i 

.i 

o,ss 

0,78 

0,78 

0,82 

1 ,0oo 

1,45 


Plate 8 shows a substantial difference between the curve deter- 
mined from the balance and the one determined by the wake method. 

V/e have to correct the values measured with the balance in order 
to take into account the Induced drag, and this corresponds to: 

= C; Ji I. 

In Table 8, we show the values of measured drag coefficients 
measured with the balance, after correction for the Induced drag: 


TABLE 8 


I® 

— 2 ® 

0 ® 

! 2 ® 

1 

40 i 

6 ® 

8® 1 

100 Cx,. 

1,07 

1 

i 1.03 

1 

1,203 1 

1 

1,5 

1 

1,9 ! 

i 


In plate 8, we show that the deviation found between the curve 
determined with the wake method and that obtained by weighing is con- 
siderably reduced after corrections are made. The disagreement 
between the two tests is in part caused by the fact that the model 
was connected to the balance with a number of wires which results 
in a turbulent flow over a large part of the wing. The measurements 
of momentum loss were performed in the median section of the wing, a 
location where the flow is not perturbed. It is, therefore, normal 
that the drag coefficient determined by this latter method v/ould be 
less than the one found with the balance (Figure 17). 
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forward attachment points 



Also, Mr. Carafoli showed that the edge perturbations come from 
the lateral play which exists between the wing and the panels. This 
produces an Increase In the drag coefficient which Is a function of 
the play. 

In plate 9 , we show the polar relative to Cx determined by the 

wake method associated vjlth values of C corresponding to an Infinite 

z 

experimental aspect ratio. We also show those found from balance 
measurements. Finally, we show the polar from balance measurements 
which were corrected for the "boundary layer effect". 

The moment coefficient of the wing was determined from balance 
measurements and also from pressure measurements performed In the 
median section of the wing. The results of these measurements are 
given In Tables 9 and 10: 


TABLE 9 


1 

1° 

— 3,6 

! 

0®,50 1 

2°.50 

4®,50 

6®, 60 

8®, 60 

10°.60 j 

100 Cm, ^ 

i 

5 

10 I 

i 

15 

19,6 

24.3 

28,7 

32.3 ! 

1 


TABLE 10 


1 

j I ® 

2 ® 

4 " 

6 ® 

8 ® 

10® 

■ 100 Cm, 

1 

13,8 

19,45 

24.4 

29 

33 
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Plate 10 shows the curves Cj^ = f(l) for these tests. 

In order to obtain the corrected = f(l) curve, from moment 
measurements with the balance, it is sufficient to shift the position 
of the reference axes as in the case of the unit lift curves. 

4 . Tests of Mr. Carafoli [4] : 

In work performed before the work of Mr. Carafoli, various authors 
tried to obtain a correct experimental infinite aspect ratio, by avoid- 
ing the short circuiting of pressures on the top side and the bottom 
side. Mr. Carafoli tried to eliminate this effect by allowing the 
edge play between the wind tunnel wall and the wing to vary. The 
extrapolation of the characteristics obtained in this way for zero 
play should then furnish the results where the edge effect is Ignored. 
Nevertheless, Mr. Carafoli found a difference in slope between the 
curves determined from the balance (extrapolated for zero play) and 
those obtained by pressure measurements in the median wing section. 

The difference noted between the reference indices relative to 

these tests is equal to 0.50° for 100 C = 100. 

z 

Using formula (15) j we determine the thickness of the boundary 
layer to the right of the wing focus. By introducing this value in 
formula (23), we can determine the vararlatlon of the angle Induced 
along the span of the wing and we therefore find an average induced 
angle . 

Tests performed by Mr. Carafoli for wings with aspect ratios 
between X = 2 and X = 8 and referred to infinite aspect ratio 
(Prandtl method) show a substantial disagreement with tests performed 
between panels. • ■ - 

} 

In plate 11, we show tests and also the polar measured for plane 
flow and corrected for the "boundary layer effect". 



5 . Tests of Mr. Glrerd [ 53 » 


Mr. Glrerd tried to obtain a correct experimental Infinite 
aspect ratio by using guard tips having the same profile as the 
portion measured and which follows the incidence variations of it. 
The tips were selected sufficiently long in order that the perturba- 
tions Induced by the transverse velocity gradient in the boundary 
layer could be avoided and could not be felt at the wing. 


The angle Induced by the "boundary layer vortices" was calcul- 
ated for a certain number of points along the span of the wing. The 
variation of this angle Induced along the measured part is shown in 
plate 12. An average Induced angle was derived from this curve. By 
assimilating the wing into a wing with finite aspect ratio it is 

possible to calculate the equivalent fictitious aspect ratio obtained 
in the test of Mr. Glrerd: 


Conclusions : 


for 


C. = 


— ‘ i i = -^ 57,3 

Ai„, for the measured sec. 

57,3 


'kf = 


« X 0,03 


= 238. 


100 , 

= 0»,08, 


The calculation method used may be Imperfect but allows one to 
obtain corrections which agree well with experiments. 

The results obtained are practically Independent of the geometri 
aspect ratio of the wings used. Under these conditions, it can be 
stated that the aerodynamic characteristics determined after the 
"boundary layer" corrections have been made are those which corres- 
pond to an infinite experimental aspect ratio. 



Chapter III 


THEORY OP A WING WITH FINITE SPAN 
I . Summary of the Prandtl theory 

Let us consider flow around a wing with finite aspect ratio. 

Aerodynamic phenomena to which the wing is subjected are not 
studied for the period In which the flow regime at velocity Is 
established. The motion of the fluid is assumed to be permanent. 

The simplified hypotheses of the linear Prandtl theory are the 
following: 

- An aircraft wing can be assimilated Into a system of vortices 
making up the v;all of the wing and connected to It and another sys- 
tem of vortices, "free vortices", which make up an infinite sheet 
(a sheet which develops in the wake of the wing). The bound vortices 
are parallel to the general direction of the span. For wings' for 
which the span is large with respect to the depth, it is assumed that 
the velocity induced by the free vortex sheet is constant along the 
depths of the wing at a given section. The elementary force of this 
section is proportional to the sum of the elementary circulations dis- 
tributed over the section. 

Therefore, one can concentrate all of the circulations at their 
center of gravity. The geometric location of these points is a single 
line which can replace the wing. The lifting surface of the wing is 
therefore replaced by a "lifting line". 

The free vortices detach at the trailing edge perpendicular to 
the span. These vortex strips essentially follow the direction of 
the general flow. A flat vortex sheet extends from the lifting line 
to infinity downstream. 



The vortex Intensity of a strip of width dy is equal to the 
variation of the circulation along this lifting line at this point: 


The free vortices induce velocity distributed along the span. 

These are induced velocities which modify the position of the refer- 
ence axes (position which was tied to the direction of the air flow 
at infinity upstream) . 
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The change in the direction of the reference axis results in a 
component of resistance (Induced drag) with the following value: 

The aerodynamic global coefficients of the wing are: 


^z=?yor^rdy, 

•/ — 

2 

/•+2- 

R'l = ? I ^^wdy — Rx,.. 

2 

II . Calculation of induced velocities (Glauert method [6] : 

Let us consider a vortex strip of intensity F which escapes at 

y 

the trailing edge of the wing at a distance y from the median section 
(Figure 18) : . 




The velocity Induced at the point M located at a distance of 
from the median section and for a strip which extends to Infinity 
downstream Is : 


dwi = — 


ix 

2 2 n 


(y — Ui) 


1 er dy 

■i-By(y — yi)' 


The velocity Induced at the point M for all strips emanating 
from the wing Is: 






cF dtf 

sy{y — yi) 
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The angle Induced In the section under consideration Is: 

, L 

't'l _ 1 dij 


yi) 


We are dealing with an Integral-differential type of system. 
Glauert gave a solution. 

The circulation over the wing is expressed In the form of a 
trigonometric series by setting: 


L 

y = — ^ ccsv 


and 


dy = ^ sin^d'^. 


L Is the span of the wing, ij; is a variable whose limits are 0 
and 7T for values of y equal to and 


from which 


and 


r =,2 L\'o y An sinn-;. 


dT = 2 LVo V n A„ cos n d 


( 1 ) 


w 


= -L 

ySyy — y, 


y nA„ cos n d 

= / J 

/ L, . 

t,/ 0 2 




y A„ cos n d 


H3 


0 


COS ^ — cos V, 



One 


The resolution of the integral (1) was given by Glauert . 
obtains 


( 2 ) 



-ri- y n\n sinn -^. 

sm — - 
1 


The theoretical circulation at a point of a wing with infinite 
aspect ratio is : 

(3) r = 4naVoSini„^7t/(l + r,e) Voiat 


a, radius of the generating circle; 

1 j aerodynamic incidence of the wing; 

1, chord of the profile; 

Hj .coefficient which varies between 0.6 and O.J3i 
e, relative thickness of profile. 

/3^ 

The aerodynamic incidence at a given section of the wing with 
finite span, is 
{^) i^ = ij — 9^ ig — 

’ 0 

where is the angle Induced at the section under consideration. 


By combining equations (3) and (4), we have: 


(5) r;r/(l -r oe) — 

V . / <r 4 YjV, ay y — y, 


In (5) v.-e replace the Integral by its value (1) 


where 


An — 


The lift component has the following expression: 


R, = pV„*L- 


f L- / ^ sin n'i» sin d ’{» 

Jo 


= ^Vo*nLSA.. 





The unit lift coefficient is: 



" Ai 


U 

S 


77 \iX 


n A| X i 
1 


77 fli xi ; 


for the slope of the unit lift curve we have: 


— .. OiA. 

di 


By comparing the theoretical lift curves for infinite aspect 
ratio and aspect ratio X, we have: 


Cr. 


77 rti 


. •'^1 


C:. 2:7(1 2(1 + 7) e) i 


or 


I- = 


C. 


2 77 (1 ^r.e) 


and 


C, 

u = — -T ; 

77 0./.’ 


The global induced angle has the following value 


9 j 


C= / I 




c.. 


2(l+7,e)j’ 


9 , = ^(1 


The induced resistance component can be expressed by: 

L 

Kti = ? I' ^ri«dy = ?Vo*L*^ 

2 

p p 

= p W L* ^ ^ ^ V.*77L* 2 

1 “ 1 

= ^ Vo-L* Ai y ^ = R, Ax (I + S), 

2 Ax" 


with 


1 


+ 8 = 


p 

S nA^n 
Ax* 



By using the unit drag coefficient, we find: 



C^i = -^ = C,Ai{l +S), 

I s\v 

Cx, = ^ (1 + 8). 

nX 


Applications to monoplane wings 


a) Rectangular wing 


The terms In the expansion of f have odd Indices, 
series Is highly convergent, we can write: 

r = 2 LVo 1 (A, sin -r A, sin 3 + A, sin 5 -r A, sin 7 


Since the 


We can select the following particular values: 

3n 


V, = ^ = 22»,50, 


= I = 45=. 


2 V 

corresponding to the values of 


V3 — = 67°, 50. 


= f = 90°. 


^ = 0,924, 0.707. 


0,383, 0. 
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It Is possible to determine the coefficients 
can also say that: 


^l’'^3’^5 



= ii = JL 

2 L 2 X 2^ 
K 


and 



We 


1+ S = A,» + 3 + 5 4- 7 A,^ ^ 



9o _ A| — 3 A, 4* 5 Aj — / A, 
y-i yi 



In Table 1, we give the values of the coefficients defined 
above. 


For the median section, we have: 


* 
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TABLE 1 


1 

K 

V 1 . 

K, j I + ,.< 

At 

Ai 

•i! 

A: 

< 

Ki 

K 

T 

t 


[A 1 

A, 

0 

0,640 

0 

0 

0 

0 

0 

0 

0 

0 

0 

00 

0,10 

0,747 

0 314 

0,171 

0,002 

0 

0 

0.134- 

0,013 

rs^ 0 

0,165 

5,850 

0,25 

0,908 

0.785 

0,351 

0,008 

0 

rs^ 0 

0,275 

0,030 

0 

0,327 

2,850 

0,50 

1,175 

1,57 

0,541 

0,025 

0 

0 

0,425 

0,065 

0,006 

0,466 

1,850 

0.75 

1,449 

2.355 

0,660 

0,042 

0 

0 

0.518 

0,09 

t 

0,012 

0,534 

1,513 

1 

1.705 

3.14 

00 

0 

0,060 

0,009 

0,0014 

0,587 

0,10 

0 019 

0,603 

1,336 

1,5 

2,225 

4 71 

0,859 

0,090 

0,016 

0.0027 

0,675 

0,14 

0,034 

0,649 

1.165 

2 

2 750 I 6,28 

0,928 

0,115 

0,023 

0,0041 

0,729 

0,17 

0.049 

0,668 

1,077 ' 

1 

2 5 

3,260 1 7.86 

0,976 1 0,136 

0.030 i 0.0055 

0,767 

0,20 

0,063 

0,680 

1.025 1 

3 

3,780 9,43 

1,011 ‘ 0,154 

0.036 ! 0.0070 

0.794 

0,22 

0,076 

0,680 

0 990 ! 


4.290 11.00 
1 

1,038 j 0,169 
1 

0,042 ; 0,0084 
1 

0,815 j 0.24 

0,088 

0,680 

0,965 


^ = ^sin-; = 1, 

9 o = Aj — 3 A 3 + 5 A 5 — / A; 

= i (ai — 3 fl 3 + 5 Oj — 7 a,), 

7o = lA. 


The variations of these coefficients are shown in plate 13 as 
a function of v = 7 — ^ — 

Plate 14 shows the values of the coefficients A for various 
values of X. 


b) Wing v;lth elliptical clrculatior distribution 


The circulation is expressed by 


r = To 



r is the circulation for the median section of the wing with 
o 

cos \p: 


^7 



r, sin <{( = 2 L\’o Ai sin <J/, 

0. As = 0, A, = 0. 

The induced velocity has the follrvlng value: 

uf = Ai sin ■; = i = C". 

sin •;< 2 L 

Thus, the Induced velocity along the span of a wing with 
elliptical circulation distribution is constant. 


r = To v/l cos* = 
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Chapter IV 


NEW METHODS FOR CALCULATING THE 
AERODYV--IC CHARACTERISTICS OF A WING 
WITH FINITE SPAN 

(Introduction of the viscosity of the ambient fluid) 

I . General remarks 

In the theory of Prandtl, it is assumed that the span of the 
wing is large with respect to its depth. The induced velocity 
Induced by the free vortices Is constant along the depth of the wing 

Let us consider Figure 19. The velocity induced by a vortex 
element of Intensity of df and of length dl at a section M of a wing 
with a finite span varies along the depth of the wing. This velo- 
city Is defined by the electromagnetic Blot and Savart induction law 

. dr dl sin X 

dw = r — ; 

Therefore, it is a function of the distance of the vortex element 
from the point under consideration. 



Figure 20 

Concentrating the elementary circulations of attached vor- 
tices at their center of gravity amounts to Ignoring the fact that 
the velocity Induced by the vortex strip varies along the chord of 
the wing. 



We examined the results obtained by Mr. Carafoli about pressure 
measurements performed in the median section of finite span wings 
and with different aspect ratios (2 to 8). 


In Figure 20, we show the curves C = f(i) for an aspect ratio 

Zi 

X and infinite experimental aspect ratio. For the geometric inci- 
dence angle 1 and for an aspect ratio X, we have a well defined 
value of C . On the curve for the experimental infinite aspect 


ratio, the value of corresponds to an incidence angle 1^ 
can be related to i b^ means of: ^ 

O 

6i is the angle Induced in the section under consideration. 


which 


We will Introduce the term "effective incidence" and by this we 
mean the incidence angle defined by (a). 


Mr. Carafoli drew plots of pressures for each orifice of the 
pressure taps as a function of the ' effective incidence angle. All of 
the diagrams should have been identical. The deviation found between 
the various curves increases as the aspect ratio decreases and as the 
orifices under consideration approach the leading edge of the wing. 

Also, we found that for finite span wings with plane rectangular 
shapes, there is an Increase in the lift near the edge borders, a 
phenomenon which is not represented by theories derived from the theory 
of Prandtl. 

Finally, in experiments one finds two vortices which originate 
at the side edges of the wing. Certain authors assume that the sheet 
of fictitious free vortices Introduced by Prandtl could be trans- 
formed into two edge vortices, assuming a rolling up of the free 

vortices, which occurs immediately at the trailing edge of the wing. 

It ^s difficult to take into account this rolling up in calculation. 

• I 

We attempted to find replacement vortex systems for a wing with 
finite span by using the following experimental facts: 



- presence of two edge vortices In the wake of the wing; 

- edge vortices which do not have the properties of vor- 
tices In perfect fluids, and the fluid In which this occurs pro- 
duces viscous motion. 

Two systems of replacement for the finite span wing have been 
proposed. They have the following common features: The edge vor- 

tices have a circulation F' which we will define later on and they 

o 

are established In a viscous fluid. 

System 1 

The elementary circulations of the attached vortices are con- 
centrated at their center of gravity and make up a lifting line 
located essentially at the focus of the wing. 

System 2 

At a given section, the elementary circulations of the attached 
voralces are distributed over the median line of the profile which 
confirms the theory of thin profiles. 

II . Distribution of velocity In a vortex In a viscous fluid 


Together with Mr. Vlllat [7]> It is assumed that the velocity 
distribution law In a vortex In a viscous fluid Is given by 


At the Aerotechnlcal Institute at Saint-Cyr, an experimental 
configuration allows one to measure the velocities in a physical 
vortex. Two small windmills A and B (Figure 21) driven by electrical 
motors M^ and M 2 turn in the same direction and with equal velocity 
if possible. 




Velocity Is measured at a section located at equal distances 
A and B which allows one to verify the curves obtained with ohose 
determined using equation (1). These curves are given In plate 15. 

Ill . System of replacement no. 1 for a finite span wing 

1 . Hypotheses 

The circulations of the attached elementary vortices are con- 
centrated at the focus of the wing profile. The position of these 
points defines a lifting line with the same transverse dimensions 
as the real wing. The two edge vortices In a viscous fluid escape 
from the wing at the focus of the extreme profiles. They follow the 
general flow direction of velocity V and extend from the wing to 
infinity 3 and Induce velocities in the wing plane. 

A2 

2 . Determination of the circulation of the edge vortices 

In order to calcuate the velocity Induced at a section of the 
wing by edge vortices, it is necessary to determine their circulation. 

We will first start with the evaluation of a circulation in the 
first approximation for all sections of the wing. For a wing whose 
depth of chord varies or for which there is warping, the circulation 
at a point of the lifting line is defined by 
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( 2 ) 


r' 


y “ 


t, 2 


» i, di I 2 


1 j chord of the profile under consideration; 

y 

igj geometric Incidence angle of the profile referred to the 
0 lift direction; 

®, slope of the unit lift curve of the profile under study 
for experimental infinite aspect ratio. 


By using equation (2) we determine the distribution of the 
first approximation of circulation along the span of the wing. The 
maximum value achieved In this distribution defines the circulation 
of the edge vortices. 


For a wing with a constant profile and chord and which is not 
twisted (rectangular wing), the circulation of the lifting line Is 
constant along the span and is equal to 

V,/. 


p a I : ( dC, \ 


1 


s 


J 




geometric Incidence of the profile referred to the 0 lift 
direction; 

slope of the unit lift curve of the profile for infinite 
experimental aspect ratio; 
translation velocity; 
chord of the wing. 


In this case, the circulation of the edge vortices Is equal to 
the circulation which would exist at all sections of the wing if It 
had Infinite experimental aspect ratio. 

Thus, the determination of the circulation of the edge vortices 
is done by not considering the velocities which they. Induce on the 
wing. 
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3 . Calculation of Induced velocities 


A3 


The velocity Induced by a vortex In a viscous fluid at the point 
M located at a distance from this vortex Is equal to the follow!- 
according to (1): 

(3) 




2kti 


In the case of the wing, we are dealing with a seml-vortex. By 
symmetry, expression 3 has to be divided by 2; we then find: 


(4) 


4 Ttl 'l 


The resulting Induced velocity 
at M (Figure 22) Is equal to the sum 
of the Induced velocities at this 
point by vortices 1 and 2. 



Figure 22 


The velocity induced by vortex 
2 located at a distance from M, Is equal to: 

(5) 

SO that the resulting Induced velocity at M has the following value: 

( 6 ) 




The velocity diagram at the point M Is shown in Figure 23. The 
aerodynamic incidence of the local profile Is modified by an Induced 
angle, whose value at the ordinate point y Is equal to: 



5U 


Figure 23 



In particular, for the median section of the wing, the induced 
angle has the following value: 


( 8 ) 


. ('«f - 


d ff/«0 = 


4 It 


1 — e ! &•< 
L 
o 
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For the median section, we can see that the viscosity terra can 
be Ignored, x = p becomes: 


s iy^o = 


/. dCz A 
2-L 


In the case of rectangular wings X = we have: 


5 t%=o — 5/, 3 


. dC. 
2:tX ’ 


( 9 ) 


5 ;0 - A / ^ ^ A 


v;hl:h Is an expression that has an identical form as the one found 
by the Glauert theory. ' 


Experimental verifications 


a) Median section 

First of all, we should mention the test performed at the Aero- 
Technical Institute of Salnt-Cyr by Mr. Carafoll. These were done 
before pressure measurements over median sections of wings having 
variable aspect ratio (2 to 8). Also, he carried out measurements 
In the mf^dlan section of a wing placed between panels. For these 
latter conditions , the correction formulas for the "boundary layer 
eff:-ot" established In Chapter II allow one to show that the equiva- 
lent fictitious aspect ratio obtained Is practically equal with the 
Infinite experimental aspect ratio. 

If we have the curves C = f(l) derived from pressure measure- 
ments In the median section of wings with aspect ratios (2 to 8) 
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Figure 21 


and the curve C = f(i) for the infinite exper- 
Imental aspect ratio (pressures measured in the 
median section of a wing between panels) from 
deviations between these curves, it was poss- 
ible for Mr. Carafoll to derive the corres- 
ponding induced angles. According to Figure 

24, for a given value of C one can associate 

z 


the effective incidence angle Igj.* The diff- 
erence between the geometric Incidence and the effective Incidence 
is equal to the experimental Induced angle: 

S Uz = iij — itt„- 

In Figure l4, we show the variation of the experimental and 
induced angle determined by Mr. Carafoli as a- function of aspect 
ratio X of wings tested. The induced angle is expressed in the usual 
form : 


I, 


= A„. 


In plate 1-, v;e show the curves of the variation of the various 
coefficients obtained from the measurements of Mr. Carafoll from the 
Glauert theory and our replacement system no. 1. 
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b) Distributions of lift along the span 

The documents about the experimental distributions of lift 
along the span of wings with finite aspect ratio which we used do not 
all contain the aerodynamic characteristics of profiles of tested 
wings (lift angle 0 of profiles, experimental infinite aspect ratio 
slope). Therefore, this data was necessary in order to calculate the 
characteristics of wings. We proceeded as follows: 

- In Figure 25, 1 is the unit curve of lift for the profile 
located in the median section of a rectangular wing with aspect 
ratio X. 2 shows the unit curve of lift- relative to infinite exper- 
imental aspect ratio. We can write: 
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The coefficient A used Is the one derived from the measurement 
performed by Mr. Carafoll: 

C„ = + Arx ig ^ =0. 

From (10), we derive: 


c,. = c. 

•• 

C-. = 

• A 


A, 



Am C-. — 




Thus, the knowledge of the unit lift coefficient of the profile 
of the median section of a rectangular wing with aspect ratio X 
allov;s one to determine the unit coefficient of lift of the profile 
under consideration for Infinite experimental aspect ratio. 


The circulation In the median section of a rectangular wing 
with aspect ratio X Is equal to 



The circulation of the edge vortices of the wing can be written 
as follows according to our hypotheses: 


so that 

( 11 ) 






/ 

from which 
/ 



/H6 
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Therefore, It Is not necessary to have the characteristics of 
the profiles which make up the tested wings In order to find the 
circulation of the edge vortices. 

Rectangular wings 


Wings with aspect ratio 5- 

Span: 7.62 m (reference [8]). 

In plate 16, we show the distribution of the experimental and 
theoretical lift (Glauert and system 1). 

For the calculation of the lift performed by using the system 
of replacement no. 1, the factor t In formula (7) which has the 
dimension of a time Is determined so that the curve calculated 
follows the experimental curve as well as possible. For this wing, 
t is found equal to t = 35Q0. 

Span: 0.90 m (reference [9]). 

The experimental and theoretical curves for this wing are 
shown in plate 17 with t = 31. 5- 

Remark . Comparison between the experimental lift distributions 
referred to the same geometric Incidence angle of a median profile for 
these two wings (plate l8) shows deviations between the curves near 
the edges. This therefore shows that the lift of the extremities 
of a rectangular wing seems to be related to the geometric dimensions 
of the wing. 

Wing with aspect ratio 3* 

Span: 4.56 m (reference [10]). 

In plate 19, we show the theoretical and experimental curves 
for t = 3000 . 



The agreement between the experimental curve and the one 
obtained from the system 1 seems to be poorer than for the preced- 
ing v/ings. In any case, the deviations found do not seem to be 
larger than those obtained between the Glauert curve and the exper- 
imental curve . 

Rectangular wing with rounded edges 

Wing with aspect ratio 6 (reference [11]). 

Span: 0.732 m. 

Plate 20 shows satisfactory agreement between the curve obtained 
by using system 1 and the experimental curve t = 27. 


Elliptical distribution wing 

V/lng with aspect ratio 6 (reference [11]). 

Span: 1.166 m. 

/A 7 

In plate 21 we show the theoretical and experimental curves for 
this wing t = 209. 


We have seen that for a wing with an elliptical distribution 
of lift along the span, according to the Prandtl theory, the velo- 
city induced is constant along the span and equal to 


S' / 

r — r IsJs. 


for the median section, we have: 


“L 

V, 


• dC, , 
4L ’ 


5 = 57,3 


• dC, 

- di 
4L ’ 
I 


— , dC. 

0/,3-jrX 

. _ 8 ig-i) di 

i 4 L ^ 


/ 


the coefficient A derived from our system 1 is equal to: 
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Ai = 


o dCi 

o7,3 -jf ao 
di 

6,28 y 


so that the ratio between the coefficients A,, 


and is equal t: 


^ = -i- = 0 64 

Ai 6,28 ' 


Remark . Various comparisons made between the curves obtained 
using our replacement system and experimental curves show that the 
approximation obtained by using the proposed calculation method is 
good. At the edge it gives the increase in the lift found in 
experiments for wings having a rectangular plan form. 

In Figure 22, we show the values of the semi-spans of wings 
prcccsed as a function of the logarithms of the corresponding t factors. 
We obtain a curve which has an Interpolated part, for which there are 
no tests in the region of the curve under consideration. 

IV. Replacement system no. 2 for a finite span wing 

1 . Summary of the theory of thin profiles 


For a thin profile with a small curvature, it is possible to 
replace a wing profile by its average line. Over this line the 
vortex layer is distributed. Let y be the circulation of an element 
of this layer with the abscissa x, for which the circulation equals: 

= (r circulation around profile) 


The elementary force at P (Figure 26) is equal to: 

/ ds = p V V ds ; 

V, velocity of the point under consideration. 




6o 



The distribution of the elementary circulations is related to 
the shape of the average line in the direction of the translation 
velocity V^. This distribution is determined by the following geo- 
metric condition: The velocity resulting at a point of the average 

line (velocity coming from the compositicr. of the velocities 
induced by the vortex elements distributed over the line and the 
translation velocity V) has to be tanger.-: any point of the aver- 

age line. 

In order to facilitate the calculations of velocity, 
induced by the vortex elements, the elementary vortices y are replaced 
by their projections on AC so that an element ds of the average line 
corresponds to a length element. dx of the line AC (Figure 26). 



The elementary speed Induced in m, projection of M, by the 
element located at P., has the following value: 


du = 


1 a r d.v 


2Tt(.t — g 2ra.r.i- — c* 

The total Induced velocity in m by elementary vortices distributed 
over AC is : 


^ l_ r'ar dx 

" 2nj^a.r.c--’ 


The condition for tangeney of the'- velocity vector at all points of 
the average line can be written as > 


( 12 ) 


^ r t . • 


dx 
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One usually expands the distribution of the circulation along 
the average profile line (Glauert) in a series In order to calcul- 
ate V. We then have 

— dx = Vq I I^Ao ( 1 cos 0) + V A„ sin n fl sin 9 jj dd. 


With 


X= ^(1 — COS0), 


v(x) = — Vo I^Ao + ^ A„ cosn 0 , 


^ = A, — i + V, A„ cosn 0. 
ax — ^ 

1 


The coefficients A„) are equal to 


Ao 




^ cos 0 d 0 . . . , A 
dx 


•=-r 

0 


^ cosn Orfn. 
dx 


It Is therefore possible if we have a given average line to 
determine the distribution of the elementary circulations along this 
line and In this way to find the lift of the profile. 


Also, If the circulation distribution is according to a selected 
law, the geometric condition of tangency of velocity vector at any 
point of the average line will allow one to determine the shape of 
this line. 

2 . Action on a thin profile of a vortex located outside of it 


a) Action on a circular contour of a vortex located outside 
of this contour. 

Mr. Pistolesl [12] showed that the circulation around a circular 
contour is modified by a vortex located outside of it. The comple- 
mentary circulation due to the Influence of an exterior vortex on 
the generating circle of a wing profile Is equal to the following 
product: , 

Ay = 47ra Awi, 


from which: 



a, is the radius of the generating circle, 

The velocity Induced can 
be written as: 


Ap = 


r'n 

2 TCf ’ 


r = the distance of the vortex 
frcm the point under consider- 
ation. For a = 0, we find 


2 rtf f ’ 




Figure 27 


b) Action on the average vortex line of a profile from an 
outside vortex. 


Let us simulate a vortex element y of the average line into 
a tlrcle with diameter dx (Figure 28). 

Due to the fact that there is a vortex T’ outside of the pro- 

o 

file, the circulation y of an element is modified by the following 
am runt : 
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with 


Ap = 


Av = 4 It ^ A p. 


r'o dx _ r„dx 

2rt(X-i)’ ^ 2(2n(X-.r)J X-x 



r 


Figure 28 
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The circulation of the vortex element then becomes 


Yi = Y + Ay, 
r'air ar , r„dr 


X rtOX « A , , 

Vi = Y = v’ - = + 


X — I dx X — i’ 
and this new elementary circulation can be defined by: 

Yi = -^dr. 

where is the new distribution of the circulation along the aver- 
age line of the profile. 


The velocity Induced by the vortex of intensity located at 
X and Induced at point m is equal to: 

dv' = — = ? — 

2 It (i — 5) 2 R 5 X (i — 5) 

The velocity induced at m by vortices distributed over AC is 
equal to: 

~‘d r, dr 




ax (X — 5) 


Finally, one obtains an integral-differential system. 


We will use a direct and approximate calculation procedure. 


In the replacement system no. 1 for the finite span wing, each 
profile of the wing operates independently of the adjacent profiles. 

Its circulation is defined by the resultant of the translation velo- 
city and the velocity Induced by the edge vortices. 

In analogy with this model, we assume that each vort^ex element 
of the average line of a profile will operate Independent from 
adjacent elements. 

51 

In this way, by considering Figure 29, the complementary Induced 
velocity. Induced at m, will have the following value: 

r' 

/;c i-fi . 

2 R (X - 5) 
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m 





- 

X * 


Figure 29 


The total Induced velocity at m Is equal to; 

Vc + 0 = o'. 


The conservation of the shape of the average line of the 
profile Imposes the condition that: 


— i = ^ = + y 

Vo dx Vo Vo 


i. 


Therefore, It Is possible to assimilate the action of a vortex 
outside of the average line of a profile at a point of the profile 
Into a local variation of the direction of the translation velocity 
This corresponds to a curve of the aerodynamic field near the pro- 
file. 


3 . Replacement system no. 2 of the finite span wing (hypothesis) 


The real wing Is replaced by a lifting surface consisting of 
vortex strips arranged parallel to the span of the wing. These vor- 
tices are "bound” to the wing. The circulation of an elementary vor- 
tex varies along the span for variable chord and twisted wings . For 
re jtangular untwisted wings, the circulation of the elementary vortices 
Is constant along the span. 

752 

I 

At a given section, the wing profile is replaced by Its average 
line. The vortex elements which make up the attached vortices are 
distributed over this line In conformance with the theory of thin 





Figure 30 


profiles (Figure 30 ). 

At the edges of the wing, the elementary vortices are super- 
imposed In order to create two edge vortices with circulation 
and ~hls circulation Is defined just like In the case of replace- 
ment system no. 1 . 

The edge vortices have the properties of viscous flow vortices 
and ohelr origin Is fixed at the foci of the profiles of the v;lng 
extremities . 


4 . Calculation of the velocities Induced at a point of a finite 
span wing 


The velocity Induced at a point A by a vortex element [perfect 

fluid of length OM and circulation F 

o 

(Figure 31)] Is given by the Blot and Savart 
formula : 





471/ p ' 


Figure 31 


since 


cos 9’ 


Vc 


= 11 ^ r 


r' 

cos 9 d 0 = ^ sin 9 o. 


Let us consider a rect^angular wing (Figure 32 ). Two edge -vor- 
tices ±F’^ escape from the wing at points and M2 such that: 

MoMi=M'oM, = |; 
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1 = constant cord of the wing. 
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Figure 
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The velocities Induced by vortices 1 and 2 (perfect fluid) at 
the point A whose ordinate is y with respect to the median section 
of the wing and whose abscissa is x with respect to the leading 
edge are equal to: 

- velocity Induced by vortex 1 at A: 
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(1 — sin0i): 


- velocity induced by vortex 2 at A: 


V2c 



[1 — sin 6,]. 


The resulting Induced velocity at A is equal to: 


Vc — Uic -f- V2c = 


r' 

•I* ft 




[1 — sinOi] + 


r' 

•L A 




[1 — sinOj] 


(13) 


_ r'o ( r 1 — sin 0;! , — sin 0/ 

4-1 I 


In particular, for the median section the velocity Induced has 
the following value: 

( 1 ^ ) r'n r 

Vc = [I — sin0o). 

- - ■ r- - 
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The velocity Induced at A by vortices 1 and 2 In a viscous 
fluid has the following value: 



It Is therefore possible to calculate the resulting Induced 
velocity. Induced by edge vortices In a viscous fluid, at any point 
of a section of a finite span wing. 

At one. point of a profile, the action of the velocity Induced 
by an outside vortex can be assimilated Into a local variation of 
the translation velocity direction (action on a thin profile of an 
outside vortex). It has the following value: 

/5il 


3y replacing v by Its value from (15) > we determine the local 
Induced angle v/hlch Is equal to: 


(16) 



In this way, we have defined a local Induced angle which varies 
at all points of a given section of a finite span wing. 


The 

pressure 

diagrams 



above calculation method allows one to plot diagrams of 
for sections of finite span wings. In order to plot these 
, It is necessary to have profile pressure diagrams which 
make up the wing, established for experimental 
Infinite aspect ratio and for different Incidence 
angles. In effect. It Is then possible to plot 
C = f(l) curves (local pressure coefficient) for 
Infinite experimental aspect ratio. Let (Figure 
33) us consider the curve obtained in this way for 
a point with abscissa x. When the profile 
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corresponds to a wing with aspect ratio X at the point under con- 
sideration, there exists an Induced angle 61 t_„ to which corres- 
ponds a pressure coefficient variation AC^ such that: 


ACp ^ 

811-1 di ’ 


dCp. 

ACp = — :j===-'SiW. 
di 


The value of the pressure coefficient at the point being studied 
and for aspect ratio X Is therefore equal to: 



Prom = /(t). curves, one can determine the characteris- 

tics Cp^ = /(0 (aspect ratio X) which allows one to find pressures 
at any point of finite span wings. 


5 . Experimental verifications 

Mr. Carafoli measured pressures in the median section of wings 
with various aspect ratios which varied between 2 and 8. He also 
made measurements in the median section of a wing placed between 
panels. Among these experiments, we only used the results of the 
tests with the aspect ratio two wing and the wing between panels. 

We believe that the selection of aspect ratio 2 was justified because 
of the fact that the curvature of the aerodynamic field increases 
when the aspect ratio decreases. We made the selection to verify 
our replacement system. 

We plotted experimental curves C = f(l) for various measure- 

P 

ment points (aspect ratios 2 and x = oo„), 
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In the following table we show data on the sides of the profiles 
used by Mr. Carafoli and shown In Figure 3^* The pressure tap ori- 
fices were located at the points which serve to construct the wing 
profile (sides given In percentage of the chord): 
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ff.V 


Figure 3^ 
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\y 
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5,96 

6,15 
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6,65 

6,27 
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5,77 
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6,52 

6,82 
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In formula (16), the calculation shows that for the median 
section of the wings tested, the term which takes into account the 
fluid viscosity can be ignored. 


We calculated the local Induced angles 
(median section) along the chord of a wing 
with aspect ratio 2 for various Incidences. 

We then related the C = f(i) experimental 

P^L 

curves for aspect ratio 2 with those for 
the Infinite experimental ratio curves. 

For this purpose, let us consider Figure 35. 

f(l) (X = 2) curve allows one to determine 61^ for aspect ratio 
X by calculating the angle induced at the point under consideration 
for a given geometric Incidence. Therefore, we find an aerodynamic 
operational Incidence 1' at the point under study equal to: 

l' = ig 5 I/. 



Knowledge of the C = 

P *\ 
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Using the Glauert method, we calculated the Induced angles for 
the median profile of a wing with an aspect ratio of X = 2 and for 
various geometric Incidences. 

In plates 23 to 28, we show the Cp = f(l) curves for X = 2 and 
X = 00 as well as the curves which one obtains by reducing thc~~. 

Q X 

for X = 2 to X = either using the Gauert method or by using our 

system no. 2. 

In his thesis, Mr. Carafoll also gives values of Induced angles 

corresponding to his measurements (determined from C = f(l) curves 

for median sections of tested wings). By using these experimental 

Induced angles and by using Figure 35, we plotted C = f(i - 61 -i ) 

P s 

which are shown In plates 25 and 26. 
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In plate 29, we show pressure diagrams for the median wing sec- 
tions. One wing Is between panels (infinite experltiental aspect 
ratio) and the other Is for the aspect ratio 2, reduced to X = “ 

w A 

using the Glauert method and using our system 2. The agreement 
obtained between the experimental curve and the one derived from 
system 2 Is satisfactory. 

V. Use of the replacement system no. 2 for the study of the 

Influence of lateral bounda r y layer s on t e sts between panels 

1 . Hypotheses 

The hypotheses Introduced In Chapter 11 for finding the influence 
of lateral boundary layers on tests between panels allow one to 
determine the unit lift coefficients of the profiles under study. 

For the same effective Incidence (median profile) and for a wing with 
a large geometric aspect ratio ano for a wing with a reduced geo- 
metric aspect ratio in a plane flow, substantial differences were 
found between the various pressure diagrams measured under these 
experimental conditions. It seemed to us to be possible to adapt the 
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replacement system no. 2 to plane flow tests In order to reduce the 
deviations among the various diagrams. 

The wing placed In a plane flow traverses the wind tunnel with- 
out a lateral play. It has a constant profile and chord and Is not 
twisted. 



Figure 36 


This wing Is replaced by elementary vortices with constant cir- 
culation arranged parallel to the span and having the same transverse 
dimensions as the real wing. 
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In a section the base profile Is replaced by Its average line. 

The elementary vortices are distributed over this line (real thin 
profiles ) . 

At the edges of the wing, the elementary vertices superimpose 
and generate two vortices In a viscous fluid which come from the wing 
at the foci of the profile located at the displacement thickness of 
the lateral boundary layers which develop along the walls of the wind 
tunnel. These extend to the foci considered to be at Infinity down- 
stream (Figure 36) (viscosity Is only Involved near the vortex axes). 
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The circulation of the edge vortices is defined in systems 1 
and 2 . 

The flow is produced between the walls and it is nere.oary to 
take into account the stresses which affect the fluid caused by the 
walls. These stresses are represented by image vortices. Figure 36 
gives the diagram of the equivalent vortex system . 

This is an infinite row of semi-vortices with alternate orient- 
ations with circulation and step L and with the following 

origin abscissa: 

2 . Calculation of induced velocity 


The velocity induced at a ordinate point Y by two vortices in 
perfect fluid with circulation ±T' of step L and located at a dis- 
tance y from the wind tunnel axis is equal to: 


with 


in particular, 


u 

ri . 


2 TZ u 
in — ^ 


Sin 


2L 2-11 2-Y’ 

cos — ^ — cos — 


9 = 5—8', 


sin ^ 


2:t 




cos 


— cos 


for the median section, we have; 


2rY 

L 



sin 




cos — I 


r' 

2 L 







1 


Let us consider Figure 37. The vortices are located at a 

distance of 5* from the wall which is small in most cases. The 
velocity induced in the median section is practically due to the 
vortices themselves and their first image. It is possible to locate 
two fictitious vortices at PP' along the extremities of the wing and 
whose circulation Is defined so that the velocity Induced at A 
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(focus of the profile of the median section) by these vortices Is 
identical with the one induced at the same point by the vortices 
and their Images: 






rtL 2L 


""r (!-*•) 

I’d-**) 


. 2x/L 
sm 


r' = - r' 


(2 -*•)-' 


cos 





We can therefore replace the wing between panels by a finite 
spar, wing (L) and one can calculate the velocity Induced by the edge 

I 

vortices of circulation at any point of the median section using 

replacement system no. 2. 
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In particular, the velocity induced at A' (Figure 37) is equal 


to : 


Uc,. 


IlsL. 

4:tL 


(1 — sin %). 


and by replacing 1" by its 


value 

2t./L 


Vc^- = 


r' 

2L 


sin 


cos 




(1 — sine,)): 


In this way one can calculate the Induced angle at any point of the 
median section of a wing between panels. At the point A’, i~ has 
the following value: 


(17: 




Sl/-A' = 



2LV, 




(1 — sin0o) 


(1 — sin 0o) 
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3 . Experimental verifications 

We performed pressure measurements in median sections of wings 
in a plane flow with geometric aspect ratios of X = 0.55 ani X = 

o o 

3.6 (tests described in Chapter II). The tests performed wi~h X = 

o 

8.6 can be considered as corresponding to a wing with Infinite exper- 
imental aspect ratio. We plotted the curve C = f(l) for aspect 

P 

ratios X = 8.6 and X = 0.55 for various pressure tap orifices 

s s 

located on the top side and the bottom side of the profile. 
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We calculated the angle Induced at these points by using system 

no. 2 of replacement for the finite span wing adapted to the plane 

flow tests for various values of the geometric Incidence angle of the 

wing. We then plotted Cp^ = f(lg-61^) curves for the aspect ratio 

X = 0 . 56 . These curves are very close to the ones obtained for the 
S 

Infinite experimental aspect ratio. 


The sides of the measurement points are represented In the fol- 
lowing table: 
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1 
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11 

12 

13 

14 

15 

16 

17 
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19 

X 

82,4 mmi 80 
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35,7 

23.8 ' 

14 i 

1 

4,5 j 

y 

3,7 mm, 0 

0 

0. ; 

0 ! 

0 

0 

0,5 ! 

1 : 


The curves C = f(l) are shown In plates 30 to 3^ 
P 
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For a given value of the geometric Incidence for a wing with 

asoect ratio [X = 0.55 (0 incidence)] (Figure 38 ) and for a given 

value of the coefficient C , the deviation betweeen 

the experimental C = f(l) curve for wings having 

P 

aspect ratios X^ = 8.6 and Xg = 0.55 represents the 
value of the local experimental Induced angle at 
the measurement point considered. 



In plate 35, we show the variation of the 
experimental and theoretical local Induced angles [derived from for- 
mula ( 17 )] along the wing chord. 


The unit lift curves of the profile determined from pressure 
measurements In the median section of wings having aspect ratios of 



X =8.6 and X = 
g g 

dence angle 1^=0 

to be equal to 2 °. 


0.55 are shown In plate 4. For the geometric Incl- 
, the effective experimental Induced angle Is found 
We show this value In plate 35. 


Finally, In plate 36 we show the pressure diagrams found In the 

median section of wings with aspect ratios of X = 8.6 and 0.55 

§ 

(geometric incidence) as well as the diagram obtained from the one 
for an aspect ratio of X = 0.55 and reduced to the infinite experl- 

o 

menoal aspect ratio using system no. 2. In this plate we also show 

the pressure diagram measured for the aspect ratio X =0.55 and 

g 

reduced to an infinite experimental aspect ratio, taking into account 
the effective Induced angle for the section under consideration 
[angle derived from the C„ = f(l) curves] 

The agreement between the experimental diagram measured for 

X, = 8,6, X, = 5c„, and the one obtained from X =0.55 reduced to the 

g 

infinite experimental aspect ratio and by using system no. 2 seems 
to ce good. 


Con e lusions 


The two finite span wing replacement systems proposed seem to 
give good results compared with experiments, as far as the distribu- 
tion of lift along the span and the distribution of the pressure 
along a given section of the wing are concerned . (for finite span or 
betv:een panels). Therefore, it is possible to follow, more precisely 
than using the Prandtl model, the aerodynamic local operation of a 
wing. This will allow one to develop a simple replacement model for 
a propeller. 

In more general terms, the Prandtl replacement system agrees 
with the Helmoltz law about vortices and disagrees with experience, 
because of the existence of a free vortex sheet. It is difficult to 
imagine the rolling up of the edge vortices. Our replacement sys- 
tem,3 on the other hand seem to agree well with experimental data 
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and disagrees with the Helmoltz laws. 


On this topic. It should be noted that the edge vortices of a 
finite span v." -g violate the fourth theorem of Helmoltz about the 
constant intensity of the vortices. 

/6l 

A more detailed study on the formation of the edge vortex of 
the wing as well as Its deterioration due to viscosity would be 
necessary to derive more exact laws than the ones discussed in this 
chapter. In particular, the correction term which involves the 
ambient fluid viscosity which we have adopted certainly will not 
only be the function of the geometric wing dimensions. 



I 
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